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FULL DISCRETISATION OF SEMI-LINEAR STOCHASTIC WAVE EQUATIONS 
DRIVEN BY MULTIPLICATIVE NOISE 

RIKARD ANTON*, DAVID COHEN', STIG LARSSON* * , AND XIAOJIE WANG § 


Abstract. A fully discrete approximation of the semi-linear stochastic wave equation driven by multiplicative 
noise is presented. A standard linear finite element approximation is used in space and a stochastic trigonometric 
method for the temporal approximation. This explicit time integrator allows for mean-square error bounds indepen¬ 
dent of the space discretisation and thus do not suffer from a step size restriction as in the often used Stormer-Verlet- 
leap-frog scheme. Furthermore, it satisfies an almost trace formula (i. e., a linear drift of the expected value of the 
energy of the problem). Numerical experiments are presented and confirm the theoretical results. 

Key words. Semi-linear stochastic wave equation. Multiplicative noise, Strong convergence, Trace formula. 
Stochastic trigonometric methods, Geometric numerical integration 
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1. Introduction. We consider the numerical discretisation of semi-linear stochastic 
wave equations of the form 

dii — Audt = f(u)dt + g(u)dW in 0x(O,°°), 
u = 0 in dS> x (0,°°), (1.1) 

w(-,0) = MO; m(-, 0) = vo in 

where u = u{x,t ) and @ C R rf , d = 1,2,3, is abounded convex domain with polygonal bound¬ 
ary dS>. The denotes the time derivative 4. Assumptions on the smoothness of the 
nonlinearities / and g will be given below. The stochastic process {W(f)} ; >o is an Li{Si)- 
valued (possibly cylindrical) Q -Wiener process with respect to a normal filtration {&t}t >o on 
a filtered probability space The initial data uq and i'o are i^o-measurable 

random variables. We will numerically solve this problem with a linear finite element method 
in space and a stochastic trigonometric method in time. 

We refer to the introductions of ED and a for the relevant literature on the spatial, 
respectively temporal, discretisation of stochastic (linear) wave equations. Further, the recent 
publication l22l presents a full discretisation of the wave equation with additive noise: a 
spectral Galerkin approximation is used in space and an adapted stochastic trigonometric 
method, using linear functionals of the noise as in El, is employed in time. Furthermore, 
the time discretisation of nonlinear stochastic wave equations by stochastic trigonometric 
methods is analysed in EQ. Finally, let us mention the recent publication El which analyses 
convergence in L P (Q .) of the stochastic trigonometric method applied to the one-dimensional 
nonlinear stochastic wave equation. 

In the present publication, we prove mean-square convergence for the full discretisation 
to the exact solution to the nonlinear problem (11.11) . Furthermore, using this result, we derive 
a geometric property of our numerical integrator, namely a trace formula. The trace formula 
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(the linear drift of the expected value of the energy) for the exact solution of ( ll.lt as well as 
for the finite element solution and the completely discrete solution are presented. 

Strong approximations of stochastic wave equations are relevant in many real applica¬ 
tions. For example, let us consider the motion of a strand of DNA floating in a liquid as 
presented in J9) and references therein. The motion of the DNA molecule may be modeled 
by a wave equation and the impact of the fluid’s molecules may be modeled by a stochastic 
force acting on the string. When two normally distant parts of the DNA get close enough, 
biological events, such as release of enzymes, occur. It is thus of interest to consider strong 
approximation of stochastic wave equations in such a situation. 

The paper is organised as follows. We introduce some notations and mention some useful 
results in the next section. Section [3] presents a mean-square convergence analysis for our 
numerical discretisation. A trace formula for the exact and numerical solutions is given in 
Section 01 Finally, numerical experiments illustrating the rates of convergence and the trace 
formula of the numerical solution are given in the final section. 

2. Notations and useful results. Let U and H be separable Hilbert spaces with norms 
||-||f/ and ||-||tf respectively. We denote the space of bounded linear operators from U to H by 
,H), and we let //MU,Ii) be the set of Hilbert-Schmidt operators with norm 

\ 1/2 

\\Te k \\l\ , 

where {e,t}“ =1 is an arbitrary orthonormal basis of U. If H = U, then we write jSf(f/) = 
2z f(U,U) and HS = .2Y(U ,U). Let Q £ J >?(!/) be a self-adjoint, positive semidefinite opera¬ 
tor. We denote the space of Hilbert-Schmidt operators from Q l ^ 2 (U) to H = U by 2z? 2 ° with 
norm 


\\ T \\j? 2 (U,H) ■= ( Y 


U-=l 


llrll^o = ||t<2 1/2 ||hs. 

For the stochastic wave equation ( Q >, we define!/ and denote the L 2 (^)-normby 

ll'll |H|l 2 (®)- Further. we set A = —A with D( A) = H 2 {<2>) nH<j(^). 

Let (£2,^,P,{J^r};>o) be a filtered probability space and L 2 (£l,H) the space of H- 
valued square integrable random variables with norm 




Next, we define the space H a = D(A n/2 ), for a £ K. with norm 

IMIa := ||A“ /2 v||l 2 (0) = (l A7(v, - 

where {(A,. ( Pj)YJ \ are the eigenpairs of A with orthonormal eigenvectors. We also introduce 
the space 


H a :=H a xH a ~\ 

with norm |||v|||^ := ||vi||^ + ||v 2 ||^_ 1 , for a € R and v = [vi,V 2 ] r . Note that H° = U := 
L 2 (&) and H := H° = H° x H . In the following we denote the scalar product by (•,•) = 
(■,-)l 2 (@) anc * recall the notation for the norm ||-|| = ||-||l 2 (@). 
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Denoting the velocity of the solution to our stochastic partial differential equation by 
«2 := M] := u, one can rewrite ( 11. Il l as 


dX(t)=AX(t)dt+F(X(t))dt+G(X{t))dW(t), t > 0, 
X(0)=X 0 , 


where X := 


U\ 

u 2 

,A:= 

' 0 I 
—A 0 

, F(X) := 

0 

/(«l) 

.G(X):= 

0 

_g( u l). 

and Xq := 

1-1 

O O 

S > 

i_i 


. The 


operator A with D(A) = // 1 = // x H° is the generator of a strongly continuous semigroup 
of bounded linear operators E(t) = e'' 1 on H = H° = H° x H 1 , in fact, a unitary group. 

Let {5j}bea quasi-uniform family of triangulations of the convex polygonal domain S> 
with /zk = diam(/T) and h = max^^ /zk- Let V), C H { \(J2>) = H l be the space of piecewise 
linear continuous functions with respect to .% which are zero on the boundary of and 
let P/h : H° —> Vi, denote the //°-orthogonal projector and S&h : // 1 —> V), the //^orthogonal 
projector (Ritz projector). Thus, 


{&hV,w h ) = (v,w h ), (V& h u,Vw h ) = (Vu,Vw h ), Vv GH°,uG H\w h G V h . 


The discrete Laplace operator A/,: V* —> Vh is then defined by 
(A h v h ,w h ) = (Vv*, Vw A ) Vw/, G V ft . 

We note that /5?/, = A /; 1 /^A- We also define discrete variants of ||-|| a and //“ by 

\\vh\\h,a = ||A“ /2 v*||, \’h e V/, 

and //“ = V h equipped with the norm || j|/,, a . Finally, the finite element approximation of 
CD can then be written as 


du /,,1 (t) + A h u h .x ( t)dt = & h f(uh,\{t))dt + P? h g(ui hl (t)) AW (/), t > 0, 
M/,,1 (0) = Uhfil M/,,2 (0) = V/, i0 , 

or in the abstract form 

d X h (t) = A h X h {t)At + 0> h F{X h {t))At + &> h G{X h {t))AW{t), t > 0, 
X h (0)=X hfi , 


where A/, := 


0 

—A/, 


X/, := 


M/,,1 

M/,,2 


F and G are as before, and A/,o := 


M/,,0 

V/ 7 ,0 


with 


m/,,0 = ^/,m o, v /,,0 = ^/,V0 G V/,. Note the abuse of notation for the projection £?hF{Xh) = 
(0, • f X > hf(uh,i )) T and similarly for ,^G(A/,). This will be used throughout the paper. Again, 


A/, is the generator of a Co-semigroup Eh(t) = e tAh on Hi, := Hj} 
We study the equations (12. U and (12.3b in their mild form 


xH7 


X{t)=E(t)X o+ f E(t-s)F(X(s))As + [‘E(t -s)G(X(s))AW{s), 

Jo Jo 

X h (t)=E h (t)X, h 0 + f E h (t-s)0> h F(X h {s))As + f E h (t-s)0> h G{X h (s))AW(s), 
Jo Jo 

where the semigroups can be expressed as 


m 


E h (t) 


C(t) A -G 2 s(t) 
-A l / 2 S{t) C(t) \ ’ 

C h {t) A h l/2 Sh(t) 

-A l „ /2 S h (t) C h {t) 


(2.4) 

(2.5) 


( 2 . 6 ) 

(2.7) 
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with C(t) = cos(fA 1//2 ), S(t) = sin(fA 1 / 2 ), C/,(f) = cos(fA / 1 /") and S/ft) = sin(fA*/ 2 ). 

In order to ensure existence and uniqueness of problem ( 11.11) we shall assume that uq £ 
LiJQ. If 1 ') and vo £ L 2 ifl.Hl *), with 7 = max(/3.1) for some regularity parameter /j > 0, 
and that the functions /: L2{3>) -G L 2 {S>) and g: L 2 {S>) —> Jzf 2 ° satisfy 

ll/(“)-/(v , )ll + lk(“)-^(' , )ll<C|| m-v||, if P >0, 

Wu)\\ + \\AV- l V 2 g( u )\\^<C(l + \\u\\), if 0</3< 1, ( 2 . 8 ) 

||A^- 1 )/ 2 /( m )|| + ||A^- 1 )/ 2 g( M )||^ < C(1 + IIA^-D/^H), if /3 > 1, 

for all m, v £ Loft?) in the first two inequalities and for all m £ H& ~ 1 in the last one. Through 
the text, C (or C\ ,Ci,K\ etc.) denotes a generic positive constant that may vary from line 
to line. We assume that the order of initial regularity 7 > 1 so that the discrete initial value 
m/, o = &hu 0 is well defined. 

Lemma 2.1. Assume that uq £ L 2 (£2,//!'), vo € with 7 = max(j3,1) and 

the functions f and g satisfy (12.8b for some p > 0. Then there exists a unique solution to the 
stochastic wave equation (ED and the finite element equation (El given by the solution of 
their respective mild equation, i. e., equations ED and ED- 

The proof of this lemma follows from H Theorem 7.4],see also the proof of Theorem 2.1 
in l 2 ll . 

We now collect some results that we will use later on. Sketches of the proofs of these 
results are collected in the appendix at the end of this paper. 

• The error estimates for the cosine and sine operators (Corollary 4.2 in [|T 6 i ): Denote Xq = 
[mo, vo ] 7 and let 

%(t)X Q = (C h {t)£ h -C(t))uo + (A h 1/2 S h (t)^ h -A-^ 2 S(t))v 0 , 

%(t)X Q = -(Al /2 S h (t)l% h - A l / 2 S(t))u 0 + (C h {t)& h -C{t))v 0 . 


Then we have 


||^(fAo||<C-(l+f)-/J r “ 1 |||A 0 ||| r , 

HWo\\<c-(i+t)-h&-v\\\x 0 \\\ y , 


t> 0 , 76 [ 1 , 3 ], 
t> 0 , 7 £ [ 1 , 4 ]. 


( 2 . 9 ) 


These will be used to estimate the error contributions from the initial values. In order to deal 
with the convolution terms in (12.5b we single out the following error estimates. Let 

■%h(t)v 0 = (A/7 S 1 , (t) £?h — A~ 1 ( 2 S (f)) v 0 , 

<%h(t)v 0 = (C h (t)^> h - C(f))v 0 . 

Then we have 

||^(Ovo||<C-(l+0-/»^||vo||p_i, t>0, p € [0,3], 

||^(Ovo||<C-(l+f)-/t^- 1) ||v 0 || |3 _i, t>0, € [1,4]. 

• The temporal Holder continuity of the sine and cosine operators, see (4.1) in Q: 

\\(S h {t) - S h (s))A h P/2 \\j? (u) <C-\t-sf , P £ [0,1], 

II (C*(f) - C/,(i))A7 tf - 1)/2 ||_ Sf(£/) <C-\t- s|^ _1 , p £ [1,2], 


( 2 . 10 ) 


( 2 . 11 ) 


Full discretisation of stochastic wave equations 
together with its continuous version: 

||(5(0-5(,))A-^ 2 ||^ (£/) <C-|t — 

||(c( 0 -c(0)a- ( ^- 1)/2 ||^) < c-k-/- 1 , 

• The equivalence of A/, and A, see the proof of Theorem 4.4 in fl5l : This uses an inverse 
inequality, hence our assumption about the quasi-uniformity of the mesh family. 

||A“^A-“v|| 2 <||v|| 2 , «g[-±,1], v€H°=L 2 (®). (2.13) 

• The equivalence of the discrete and continuous norm, see (2.13) in Jl]: 

c\\A. r h v h \\ < ||A r v/j|| <C||a£v a || for v h £V h and ye[-±,±], (2.14) 

Using the above estimates, one can deduce the following regularity results for the exact 
solution to our stochastic wave equation (11. Il l and for the exact solution of the finite element 
approximation (12.21) . 

PROPOSITION 2.2. Let [u\,U 2 \ T be the solution to (11.11) . where the initial values satisfy 
n 0 G Li (£2,//f), vo G L 2 {Ll,Hwith y= max(/3,1), and the functions f and g satisfy (12.81 ) 
for some 0 > 0. Then it holds that 

sup E[||ni(/)||jg + ||n 2 (f)||^_i] < C 


P G [0,1], 
P G [1,2]. 


( 2 . 12 ) 


and, for 0 < s < t < T, 

E[||„ 1 (t)-ni(,)|| 2 ]<C|t- S | 2 ^’ 1 )(E[||„o||| + ||vo||^ 1 ] 

+ su p e [ 1 + IImi('-)II^])- 

r€[0,r] 7 

The proof of this proposition is very similar to the proof of Proposition ^. 31 given below 
and is therefore omitted (see also the proofs of Proposition 3.1 and Lemma 3.3 in H21 1 ). 

The next result will be useful in Section [4] when we will deal with the trace formula of 
the numerical solution. 

PROPOSITION 2.3. Let [uh\,Uh2\ T be the solution to the finite element problem (12.2b . 
where the initial values satisfy no G Li(£2,//T), vo G L 2 (Ll,HY~ l ) with y = max(/3,1), and 
the functions f and g satisfy (12.81) for some fi G [0,2]. Then it holds that 

sup E[||n/, ,1(011^ + lk, 2 (0lli/i-i] <C 


and for 0 < s < t < T 


E[||«m( 0 — « M Mil 2 ] < C|f - s \ 2mia ^ (E[||n /; ,o||^ + Holl^-r] 

+ sup E[l + ||n /u (r)||^]), 
re[o,r] 7 


where we recall that n/, o and v/,p are the initial position and velocity to the finite element 
problem. 
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a /2 

Proof. Let us start with the first estimate of the norm of A ^ u^\(t) and consider the 
expression 

Af /2 M/,,i(0 = A^J 2 Ch {t )uhfi + ^ /2 S/, (t)vhfi 

+ Af- X)l2 S h {t-r)3P h f(u h ,\(r))&r 

+ J A}f _ 1 )/ 2 5/,(f-r)%(M/ U (r))dff(r). 

Using the fact that A/, and C/,(f) commute, the boundedness of the cosine operator, together 
with our assumptions on the initial values for the finite element problem, we get 

E[\\A P h / 2 C h {t)u hfi \\ 2 ]<C for J3 G [0,2], 


Similarly, one obtains 


E[||A^- 1 )/ 2 5,(f)v /li0 || 2 ] < C. 

To estimate the third term, we use ( 12.13b , the assumptions on / given in ( 12.8b . and the equiv¬ 
alence of the norms stated in (12.14b . First for /j g [0,1], we get 


E 


j A (P ^! 2 s h (t-r)& h f(u ht \{r))dr 


< Ci + C 2 

< 


fn\M )\\ 2 

j 0 


dr 


C 3 +C 4 J E[||i< /lil (r)||^]dr, 
because Si,(t) and A /; ^ 2 are bounded. For (3 g [1,2], we have by (12.13b 

j’ S h (t - (r)) dr " 

: rE[l + ||A^- 1 )/V 1 WH 2 ]dr<C 1 +C 2 fE[ Kl ( 

Jo Jo 


<C 


^ C 3 + C 4 


'")ll*,JS-l] dr 


J m\u hA (r)\\ 2 hJ} }dr. 


Finally, Ito’s isometry, equations (12.14b and (12.13b . and the assumptions (12.8b on g give us 


E 


J o S h (t-r)^ h A (/3 1 )/ 2 A (/3 1 )/ 2 g(M/,j(r))dW(r) < C 3 + C 4 E[||i</ lil (r)||^ j j]dr. 


All together, for [3 g [0,2], one thus obtains 


E [IKi(0llL] <Ki+K 2 J o E[\\u K 1 (r)\\ 2 hl} ]dr 

and an application of Gronwall’s lemma give the desired bound for E[||m /,.1 (Oil* p}- 

The proof for the other bound is done in the same way except for a slight difference in 
the initial values and that A^ l ^~Sh(t — r ) in the integrals is replaced by A jj* *^ 2 Chit — r). 
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We now prove a Holder regularity property of the finite element solution. We write, for 
0 <s<t<T, 

- u h: i(s) = {C h {t)-C h {s))u h , 0 +A a (S h (t) - S h {s))v h ,o 

+ J A;7 1 / 2 {Si, ( t-r)- S h (s - r)) ^ h f{u h ,\ {r))dr 

+ / A /7 I/2 ^(f-'')^ , / ! /( M / ! ,iW) dr 

Js 

+ J \ l/ 2 {Sh{t -r)- S h (s - r)) 0 > h g{u h .i (r)) dW(r) 

+ ^ A / 7 I/ 2 5/ 1 (/-r)^ , /l g(M /l j(r))dW(r). 

To estimate the first term we use (12. 1 11 1 to get 

E[|| (C*(0 - C / ,(5)) M/ ,,o|| 2 ] = E[||(C/,(f) - C h {s))A- h P ' 2 A P J 2 u hfi \\ 2 ] 
<C|f-s| 2 %[||Aj%>|| 2 ], 


for /3 e [0,1]. For j3 G (1,2] we note that \ = A^^A ^ -1 ^ 2 and that A ^ _1 ^ 2 is 

bounded in the operator norm. Using a similar argument for the second term, we get the 
following estimate for the first two terms 


E[|| {C h (t ) - C h (s))u hfl + A h 1/2 {S h (t) - 5 /! (i))v^ 0 || 2 ] 
<C|f- 5 | 2min (P- 1 )E[||«* i o||2 i/J + ||v A ,o||^_ 1 ], 

for p G [0,2]. In order to estimate the third term, we use (12.1 11 1 , the assumptions on /, and 
the equivalence of the norms given in (12. 141) . First for ft G [0,1], we obtain 


E 


J \ l/ 2 (S h (t -r)-S h (s- r)) 0 > h f{u,a {r))dr 

<C\t-s \ 2 sup E[l + ||H M (t)|| 2 ] 
te[0j] 

<C\t-s \ 2 sup E[l + || M/ , 1 (f )|| 2 ]. 
re[0,r] 


For p G [1,2] we have, using (12. 1 II) . (12. 1 3b . (12. 14b and the fact that A /( ^ *^ 2 is bounded in 
the operator norm 


E 


J q A /( 1 // 2 (5/,(f — r) — Sh{s — r))ZPhf {ui, \ {r)) dr 

<£E[||(5 /! (f-r)-5 /l ( S -r))A7 1/2 A7 (/3 - 1)/2 A( /3 - 1)/2 ^A-^- 1 )/ 2 
X A (/ 3 “ 1 ) / 2 /(M / ,, 1 (r))|| 2 ]dr 
< C|f — s | 2 sup E[l + \\u hA {t)\\h-i\ 

te[0,T] 

<C|f-s | 2 sup E[l + || M/l , 1 (f)|| 2 /3 ]. 
fe[o,r] 
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Similarly we get for the fourth term 


E 


J A h l,2 S h {t-r)& h f(u, h i{r))&r 


< C\t — s\ 2mm ^' 1 ' 1 sup E[l + ||M/ U (f)||^]. 
/e[o,r] 


To estimate terms five and six we use Ito’s isometry, (12.1 11 1. (12.131) . (12.14b and the assumptions 
on g to get, for p £ [0,1], 


E 


J \ 1/2 (Sh{t-r)-Sh(s-r))^> h g(u htl (r))dW(r) 

< J^E[\\ ( S h (t - r) - S H (s - r))A / ; / 3 / 2 A|f- 1)/2 ^ 
x g(n/, i(r))||^,o] dr 

< C\t —s\ 2 P sup E[||Af ~ 1)/2 g( M/M (f)) ||^o] 

re[0,r] 2 


2-1 


< C\t ~s\ 2 P sup E[l + || M/!il (f)||^] 
tG[o,r] 


and 


E 


J \ 1/2 Sh(t-r)&> h g(ui h i{r))dW(r) 

< f E[||5/,(f - r)A h P/2 A[^ /2 ^ h g( Ulul (r)) ||^ 0 ] dr 

J S 2 

< C\t — s\ 2 P sup E[||A}f _1)/2 g(M /!il (f))||^ 0 ] 

iG[o,r] 2 

<C\t-s\ 2fi sup E[l + || M/i , 1 (f)|| 2 ^]. 
iG[o,r] 


For p g [1,2] we again use that A^ -1 ^ 2 is bounded in the operator norm. 

Collecting the above estimates give us the statement about the regularity of the finite 
element solution. □ 

3. Mean-square convergence analysis. Recall that the exact solutions to (12.11 ) and (12.31 ) 
solve the following equations 


X(t) = E(t)X 0 + f E(t-s)F(X(s))ds + [‘E(t - s)G{X(s))dW(s), 

Jo Jo 

X h (t)=E h {t)X K 0+ f E h (t-s)& h F(X h {s))ds+ f E h {t- S )& h G{X h ( S ))dW(s), 
Jo Jo 


where Xq 


u 0 
vo 


A/,,o 


■'X’hiio 
&hv 0 


and 


m = 


C{t) A~ l / 2 S(t) 

A 1 / 2 5 (f) C(t) 


E h (t) = 


C h (t) 

-A l /2 S h (t) 


-1/2 


S h (t) 


C h (t) 


with C(t) = cos^A 1 / 2 ), S(t ) = sin(fA 1 / 2 ), C/,(f) = cos(fA*/“) and 5/,(f) = sin(fA*/ 2 ). 
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The explicit time discretisation of the finite element solution ( 12.3b of the stochastic wave 
equation using a stochastic trigonometric method with stepsize k reads 

U n+1 = E h (k) U" + E h (k) &> h F ( U" )k + E h (k) & h G( U n ) AW", 


that is, 




c h (k) \ l/2 s h {k) 


u?\ 


r A]j 2 S h {k) C H (k) 



\ 1 / 2 Sl:(k) 
Ch(k) 


&> h g(u r ;)Aw" 


A* l/2 Sh(k) 
C h {k) 


■nnm)k 


(3.1) 


where AW' 1 = W (t n +i) - W(t n ) denotes the Wiener increments. Here we thus get an approx¬ 
imation UJ « u h j(t n ) of the exact solution of our finite element problem at the discrete times 
t„ = nk. Further, a recursion gives 

n— 1 n— 1 

U" = E h ( t n )U° + £ E h 0 tn ~ tj ) 9 h F ( U j ) k+^E h (t n - tj ) & h G{W) A . 

j =o j =o 

We now look at the error between the numerical and the exact solutions U n —X(t„). We 
follow the same approach as in |[23l for parabolic problems, see also sni, and obtain 

E[||t/"-X(t„)|| 2 ]<3(E[||Err 0 || 2 ]+E[||Err d || 2 ]+E[||Err s || 2 ]), 

where we define 


Err 0 : = (E h (t n )&> h -E{t n ) )X 0 , 


Err d :=V [ 1 (E h (t n -tj)^ , h F{U ] )-E{t n -s)F{X(s))\ 
j=o J ‘j v 2 


ds 


and 


Err s := 


:= £ [ ,J+ ' (E h {t n - tj)0> h G{W) —E(t n - s)G(X(s))] dW(s). 

j=o Jt J 

We next estimate the above three terms. 

Estimate for the initial error Erro. By ( 12.9b , the first component reads 

E[|| (C h {t„)S h - C(t n ))u 0 + (A ~ l,2 S h {t n )^ h - A- 1 / 2 S(f„))v 0 || 2 ] 


< C(1 + t n ) 2 h 2 ^ (E[||« 0 ||y+ Hvollr-t])', 
for y e [1,3]. Similarly for the second component 

E[||-(Ay 2 Sft(f„)^ ft - A l / 2 S(t n ))u 0 + (Cft(f„)£»ft - C(f„))vo|| 2 ] 
<C’(l+ ? „) 2 /z^(r- 1 ) (E[||« 0 || r +||v 0 || r -i]) 2 , 


for ye [1,4]. 
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Estimate for the deterministic part, Err,]. We write the deterministic error as 

Err d =£ [ 1 (E h (t n -tj)^ l h F(U i )~E(t n -s)F(X(s))\ds 
7=0 W 

= £ [^ E h [t„ - tj)& h (F(U j ) - F(X(tj)))ds 
j=0 Jt i 

+ £ [ ,J+1 E h (t n - tj)!? h {F{X{tj)) - F(X(s)))ds 

7=0 Jt j 

+ £ [ 1 C E h{U-tj)0 e h -E(t n -t j ))F(X{s))ds 
j=o Jt J 

+ £ [ J+ (■ E(t n -tj)-E(t n -s))F(X(s))ds 
j=0 J, i 

=■ h +h+h+h, 

and estimate the second moment of each term in the above equation. For the first component 
of the first term we get the following estimate by using (12.12b and (12.8I > 

(E[||/[i,i]|| 2 ]) 1/2 < £ f J+ ' (E[\\A- h l/2 S h (t n -tj)& h (f(u{) -/(«ifo)))|| 2 ]) 1/2 d,v 
j=o J, J 

H-l , N I/O 

<c£*(E[||i//-«(f ; -)|| 2 ]) , 

f=o 

so that 

E[P[i,i]l| 2 ] < («£ (E[||t//-«( f; )|| 2 ]) 1/2 ^ <Ck n fE[\\u{-u( tj )\\ 2 ]. 

The second component is estimated in the same way 

(E[||/[i i 2 ]|| 2 ]) 1/2 <£ f i+l U[\\c h (t n - tj)^> h {f{u{) - /(ill (tj)))\\ 2 ]) 1/2 ds 
j=o J, J 

n-\ . \ 1/2 

<c£^(E[||t//- M (t 7 )|| 2 ] . 

7=0 

For the second term, using Proposition 12.21 we get 

(E[||/[2.i]ll 2 ]) 1/2 

<£ f ' + (e[||A“ 1/2 5 a (?„-7 ; )^(/(mi(7 7 ))-/(mi( 5)))|| 2 ]) 7 ds 

7=0 W V 2 

<C£ f tj+ ' (E[|| Ml (0)-«i(s)|| 2 ]) 1/2 ds 

7=0 j 

<c£ / 0+1 lo-sn n(/M) «1 s(e[||m 0 ||| + llvoll^j] + sup E[1 + ||« 1 (0llp])' ^ 

Po v H ^ te[o,T] F y 

< CA ,min(/3,l) ; 
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for P £ [0,3]. Thus 


E[||/[ 2i i]|| 2 }<Ck 2 mm W’ l \ 

_ 1 /2 

The second component /p, 2 ] has the same expression as Iy\ 2 \ except that A /? ' Sh(t n — tj) is 
replaced by C/,(f„ — tj). The same estimate holds since the cosine operator is bounded. 

The third term reads, using in (12.101) and ft £ [1,3], 

(E[||/ [3i1] || 2 ]) 1/2 < '£ f j+1 (E[||(A ft - 1/2 5,(f n -tj)& h - A- 1/2 5(t„ -f,■))/(«! Will 2 ]) ^ ^ 

^ 1 rtj -|_i i 

= L / (E[||J%(4-t ; )/(Mi(i))|| 2 ]) /_ ds 

j=0 Jt J 

<C^”f f i+l (E[||A« 5 - 1 )/ 2 /( Ml (5))|| 2 ) 1 / 2 d, 

7=0 ■'O' 

<Ch& f sup E[l + || M i(f)||^ 1 ] N | 

V6[°T] / 

<Ch^. 

For j 8 £ [0,1] we simply note that 

E[||A^- 1 )/ 2 /(« 1 (,))|| 2 ] < CE[||/(m 1 (s))|| 2 ] < c. 

The estimate for the second component is done in a similar way using now in (12. 1 01 ) 

with p £ [1,4], 


E[||/ [3i2] || 2 ]<C/^- 1 ), 


For the fourth term with j3 £ [0,3], using (12.121 ) and the assumption on the function / in 
( 12 . 8 b . we get 

(E[p[4,n I! 2 ]) 1/2 < I f i+1 (E[|| (S(t n - tj) - S(t n - ,))A-V2 /(mi (s)) ||2]) 1/2 diV 

7=0 A; 

< cl [ ti+1 (II (S(f„ - tj) - s(t„ s))A^r- 1 | 2 . E[||/( M1 (,)) II 2 ]) 1/2 d, 

7=0 ' Jt i v 2 

<c£ / (|5-0| 2 E[l + || Ml ( s )|| 2 ]) 1/ -d, 

j=0 J, j 

< Ck. 

Thus we obtain 


E[P[4,t]l| 2 ]<CA: 2 . 
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For the second component we get 

(E[||/ [4i2] || 2 ]) 1/2 < £ f i+l (E[|| (C(t n tj) - C(t n - ,))/(«! (5)) II 2 ]) 1/2 d, 
j=o J, J 

<c't f i+l (mt n - tj )-c(t n - S ))A-^ 2 \\% m 

7=0 Jt i 

x E[||A^“ 1)/2 /(ni(>))|| 2 ]) 1 dIs 

- C X [ + (l s— ol 2 ( ^ _ 1 ) ®[ 1 +ll M i( ,s )ll^-i ]) 7 ds 

j=o Jt J 

< Ck min (P~ l,1 \ 


for/3 > 1. 

Altogether we thus obtain 

E[||Err dJ || 2 ]<C-(/J+^ 2 min (^ 1 )+^£ 1 E[||l//-n 1 (t J )|| 2 ]) for j3 e [0,3], 

7=0 

n— 1 

E[||Err d . 2 || 2 ] < C-+fc 2 ^- 1 - 1 )+* £ E[||I//-«i(f ; -)|| 2 ]) for /3 € [1,4]. 

7=0 

Estimate for the stochastic part, Err s . We rewrite the stochastic part as we did for the 
deterministic part of the error: 

Err s = £ f J+> CE h (t„ - tj)0> h G{W) — E(t n - s)G(X(s))] dW(s) 
j=o Jt J 

= £ f J+> E h (t n -tj)&> h {G(U j ) - G(X(tj)))dW{s) 

j=o Jt J 

+ 1 f' +l E h (t n -tj)& h (G(X(tj)) - G(X(s)))dW( S ) 

j=0 J, j 

+ E f 1+> {E h (tn-tj)£? h -E(t n -tj))G(X(s))dW( S ) 

7=0 W 

+ E f J+l ~ tj) E(t n s))G(X(s)) dW(s) 

7=0 W 

= :/( + 72 +-^ 3 +^ 4 - 

The estimate for the first term follows by using the Ito isometry, the boundedness of 5/, 
and A h 1//2 , and the Lipschitz condition on the function g in ( 12.8b 

E[||/[ M] || 2 ] = E f i+l m\\' /2s h(t n - tj)^ h {g{u{) - g( M1 (fj))) n^ 0 ]ds 

j=0 Jt i 2 

<C*£E[||g/-« i (/ ; -)|| 2 ] 

7=0 
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for p £ [0,3]. The same estimate holds for the second component J , 2 \ with ft £ [1,4]. For 
the first component of the second term, using Proposition ^. 21 we obtain 

n— 1 rtj+i 

E [ll- / [2,i]ll 2 ] = E / E [|l A / 7 1/2 5/ I (f„-o)^fe(«i(o))-^(»iW)ll>]d, 

j=0 J, j 2 

< c E I + e oi m i( ? ./)- m i( s )II 2 ] ds 

j=0 Jt i 

<c’t \>1 -s\ 2w *W dLs( E [||uo||J + ||m|||_,] 

j=0 Jt i 

+ sup E[1 + || M 1 (0||^]) 

(€[o,r] ' 

< Ck 2min (P’ l \ 


for P £ [0,3]. Similarly, the estimate for the second component of .li reads 


E [ii^[2,2]ii 2 ] = E [ J+ m\c h (t n - t^^higiuipj)) - g(ui(s))\\ 2 ^ 0 ]ds 
j=o J, j 2 

< C^ 2 ^E[||mo ||[3 + ||vo[||_i] + sup E[1 + ||Mt(f)||p]j • 

For the second component we have ft £ [1,4], so that min(/3,1) = 1. For the first component 
of the third term we use (12. 1 Ob with ft £ (1,3] to get 

e [||-'[ 3, 11 I| 2 ]=E r +lE [||(A7 1/2 5,(4-t 7 0^-A- 1/2 5(f n -f;0k(« 1 W)ll^o]d S 

j=0 J, J 2 

A ^— 1 ft i- 4 _l 

52 / E[l|A' f - 1 >' 2 s(«,W)ll 2 ]<k 

/-« J ’I 

4B 0 4B 

<Ch^~ sup E[1 + ||Mi (f)|| »] < C/;T" 
re[ 0 T] 

by Proposition 12.21 The estimate for ft £ [0, 1] is obtained in the same way. For the second 
component, we also obtain 

E[||^[3,2]H 2 ] < sup E[1 + ||w 1 (f)[|n] < Ch^~^ 

re[o,r] 

for P £ [1,4]. Finally, for the first component of the fourth term, we get 


E [ll-/[4.i] II 2 ] = E f +1 E [ll (S(tn tj) - S(t n S ))A-P/ 2 A^y 2 g( Ul (s))\\% 0 } ds 

j=0 Jt j 2 

<CE f 1 \s-tj\ 2l3 ds sup E [1 + 1| mi (f )|| 2 ] 

j=0 Jt i ;6[0T] 

< Ck 2 P sup E[1 + || Ml (?) || 2 ], 
re[o,r] 
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for/3 £ [0,1]. For/3 > 1, wenote that A P/ 2 =A */ 2 A $ ^Z 2 and that A $ x )/ 2 is bounded 
so that we get 

E[||/[ 4 i i]|| 2 ]<Cfc 2min ( /3 ’ 1) . 

Similarly for the second component, using the regularity of the cosine operator, we obtain 

E[||^ 2 ]l| 2 ]<« 2mm(/ '— U) , 

for/3 > 1. Altogether the estimate for the stochastic error reads 

E[||Err M || 2 ]<C- (h$ + * 2min(W) +k £ E[\\U{ - Hi(f,-)|| 2 ]) for j 8 £ [0,3], 

j =o 

E[||Err s , 2 || 2 ] < C- (h^ +* £ E [||t//_ Ml (^)|| 2 ]) for p £ [1,4]. 

v f=o 2 

Collecting the estimates of the three parts of the error, we thus obtain the following 
estimate for the error in the position and velocity of the stochastic wave equation 

E[[|Cf- M i(t n )|| 2 ]<C-(^+A : 2 min ^- 1 )+A:£E[[|t//- M i( 0 )|| 2 ]), p £ [0,3], 
v l=o 2 

E[||t/ 2 ”- M 2 (M|| 2 ]<C-(/t3(/3- 1 )+A : 2 ^- 1 . 1 )+A:£ 1 E[||f//- Ml (f ; -)|| 2 ]), P £ [1,4]. 
v l=o 

Using the above error bounds and an application of the discrete Gronwall lemma proves 
the following result for the mean-square errors of the full discretisation of the semi-linear 
stochastic wave equation with a multiplicative noise. We assume that y is large enough so 
that the stochastic error dominates over the initial error. 

THEOREM 3.1. Consider the numerical discretisation of the semi-linear stochastic wave 
equation with a multiplicative noise CD on a compact time inten’al [0,7'] ,T> 0 , by a linear 
finite element method in space and the stochastic trigonometric method (ED in time. Assume 
that uq £ L 2 (f2,//7), i'o £ L 2 (f2,// 1 '~ 1 ) with y > 1 +2)3/3 and that the functions f and g 
satisfy ( 12 . 8 b for some p > 0 for the error in the position (and for some p > 1 for the error in 
the velocity). Then, for t n £ [0,7], the mean-square errors read 

\\U?-u hA (t n ) ll^o) < for p e [ 0 , 2 ], 

\\UZ-u K2 {t n )for P £ [1,2], 

Ilf/r-MiWIWo) <C- (/,^+^-D) for P £ [0,3], 
ll^ 2 -« 2 (fn)|L 2 { a,H 0 ) <C.(/t^+r i ^- 1 ' 1 )) for P £ [1,4]. 

Observe that the error estimates between the finite element solutions and the solutions 
given by the stochastic trigonometric method are proven in a similar way as above, using in 
addition Proposition l2.3l 

4. A trace formula. In this section, we will only consider the problem (II.lb with ad¬ 
ditive noise (g = 1 in ( II. lb ) and the nonlinearity flu) = —V'{u) for a smooth potential V. 
We will further consider a trace-class Q -Wiener process W, i. e., Tr (Q) = HG^IIhs < °°. In 
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this situation, the exact solution of our nonlinear stochastic wave equation satisfies a trace 
formula (see for example 12 [5) for linear stochastic wave equations), where, in analogy to 
deterministic problems, the “Hamiltonian” function is defined on //’ = H 1 x H° as 

H(X) = ^ ( (|M 2 | 2 + |Vw 1 | 2 )ck+ f V(ui)dx 
= ^ II “ 2 1| 2 + ^l|A 1/2 ;/i || 2 + J^V(ui)dx. 

In this section, we restrict our attention to additive noises, since, in this case, we obtain an 
elegant and tractable expression for the drift term in the trace formulas (see below). This is 
not the case for the case of multiplicative noise as explained in a remark after the proof of the 
next proposition. 

The trace formula for the exact solution to our stochastic wave equation is given in the 
following proposition. 

PROPOSITION 4.1. Consider the nonlinear stochastic wave equation (O with additive 
noise, that is with g = 1. Further, let f(u) = —V'{u)for a smooth potential V, let {W (f)};>o 
be a trace-class Q-Wiener process, and let the Hamiltonian H be defined as above. Then 
the exact solution, X(t) in equation (12.4b . of the nonlinear stochastic wave equation (HD , 
satisfies the trace formula 

E[ff(X(0)]=E[ff(X(0))]+^Tr(<2), t > 0. (4.1) 


Proof. Indeed, using Ito’s formula (one can apply Theorem 4.17 in [ 8 ] since X(t) is an 
Ito process and the potential V is smooth enough) for the above Hamiltonian, we obtain 

H(X(t))=H(X( 0))+ / (ff , (X(s)),GdW'(s))+ f (H'(X(s)),AX+ F(X))ds 
Jo Jo 

1 f Tr[H"(X(s))(GQ 1 / 2 )(GQ 1 ? 2 )* 

2 Jo 


+ ■ 


for all time t. Here we have G = 


1 d,v 


, since we are concerned with additive noise. The 


expected value of the second term in the above formula is seen to be zero. Using the definition 
of A and of the nonlinearity F, the integrand present in the third term reads 

(Ami,M 2 ) + (V'(mi),M2) + (m 2 ,-Ami - V'(mi)) = 0. 


Finally, using the above definition of G and the fact that the operator Q is self-adjoint, the last 
term in the above formula is seen to be equal to 

\ f , Tr(Q l,2 (Q 1/2 y)ds = t 1 -Tr(Q). 

2 Jo 2 

This shows the trace formula (14.1b for the exact solution of our problem. □ 

Remark 4.2. Similarly to the above computations, for the case of multiplicative noise, 
one would obtain 

H(X(t))=H(X( 0))+ [\H'(X(s)) 1 G(X(s))dW(s))+ [\h'(X(s)),AX+F(X)) ds 
Jo Jo 

+ ^i) ,Tr[H,,(x(s))(G(x(s))el/2)(G(x(s))el/2r]dJ 

= H(X{ 0))+ [ t (H'(X(s)),G(X(s))dW(s))+ 1 - /'Tr(g(M 1 ( i ))e 1 / 2 (g(«i( S ))e I/ 2 )*)d S . 
JO 2 Jo 
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Taking expectation thus leads to 

E[H(xm=E[H(X(0))}+ ] 1 | o r E[Tr(g( Ml ( S )) 2 (g( Ml (*)))*)] ds, 

which, for general multiplicative noise, do not give a tractable expression of the drift in the 
trace formula. 

REMARK 4.3. The trace formula is also related to the energy equation, a tool that can 
be used to analyse the existence, or nonexistence, of solutions to stochastic nonlinear wave 
equations, see for further details on this topic. 

We next observe that, for the finite element solution X/„ one has 

1 f 

H{X h ) = dK,2|| 2 + d|A / 1 / 2 M / ,. 1 || 2 + / V(u hA )dx, 

2 2 

because || Vv /z || = ||A l ! 2 v h \\ = ||A fv„|| for finite element functions v/,. This results from the 

definitions of A 1 / 2 and A }/ 2 , see Section^ Using similar arguments as in the proof of the 
above result, one can now show that the finite element solution X/,(t), defined in (12.5b . also 
possesses a trace formula. 

PROPOSITION 4.4. Let f, g and W be as in Proposition 1 4. 1\ The solution of the finite 
element approximation of problem CD, X/fit) in equation ( 12.5b . satisfies the trace formula 

E[ff(^(f))] = E[ff(X„(0))] +t^Ti{& h Q& h ), t > 0. (4.2) 


We will now prove that the full discretisation of the stochastic wave equation, that is the 
numerical solution given by (13.1b . satisfies an almost trace formula. Indeed, as seen in the 
theorem below, we get a small defect of size <£? , (T mm ( 2 (/l- 1 ), 1 )) However, due to the use of 
Gronwall’s inequality, the defect term is not uniform in time. 

THEOREM 4.5. Let f, g andW be as in Proposition \4.1\ and \4.4\ Let further the assump¬ 
tions in Theorem\3J\be fulfilled with [5 £ [1,2]. Then the stochastic trigonometric method 
CD satisfies an almost trace formula 

E[//(t/")] = E[H(U o )\+t„^Tr(0 , h Q3? l h ) + (4 . 3) 


for 0<t„<T and j3 £ [1,2]. 

Proof. The proof uses similar techniques as the ones used to prove the mean-square error 
estimates for the numerical solution in Section [3] 

To prove the almost trace formula (14.31) . we first add and subtract the expectation of the 
Hamiltonian for the finite element solution X/fit) 

E[H(U n )] = E[H(U n )-H(X h (t n ))}+E[H(X h (t n ))} 

= E[H(U n )-H (**(#„))] +E[ff(X*(0))] +t n \Tr{& h Q&h) 

using Proposition l4.4l We will next show that 

E[H(U n ) —H(X h (t„))] = i).i)) 


(4.4) 
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for £ [1,2]. Indeed, we have that 


E[H(U n )-H(X h (t n ))\=E 


\ [ {m 2 ~\u hP Stn)\ 2 )te 

2 Jg, 

\ / dA;/v;f-|AyVi(^)i 2 )d-v 

l J9 

[ (V(U?)-V(u hA (t n )))dx . 

J9> 


(4.5) 


Thus we get three terms to estimate. Using Cauchy-Schwarz inequality, the first term in the 
above equation can be estimated by (neglecting the factor 4 for ease of presentation) 

Nl^ll 2 -IIMMII 2 ]! = |E[(t/ 2 «+« /Ii2 (4),t/ 2 "—« /Ii2 (/„))]| 

< (E [||Ul + u K2 (/„)|ls,p_ 1 ]) 1/2 (E[||l£ - u^Wl^f 2 

< C (E[||A[ 1 “ /3)/2 (f/ 2 " - || 2 ]) 1/2 , 

where we have used the discrete norm, the fact that the finite element solution m /,,2 (0 is 
bounded in the mean-square sense (see Proposition l2.3l >. and the fact that the numerical solu¬ 
tion given by the stochastic trigonometric method is also bounded, i. e. 

n\\unl,p+\m\lp-i]<c<« for «=o,i,...,a-i. 

The proof of these estimates is similar to the one for the finite element solution given in 
Proposition 12.31 except that we now have a sum of integrals of length k. This causes no 
problem since we can simply use the triangle inequality for the deterministic integrals and for 
the stochastic integrals we use the property that they are independent with expected value 0. 

Using the definition of the time integrator and similar techniques as in the proof of the 
mean-square convergence, one next estimates 


A (t ^\u^-u ht2 (t n ))=Y, f J+ ' A [ l M 2 (C h (t n -tj)-C h (t n -s))<?> h dW(s) 

j=0 Jt i 

+ E /' A h l ^ / 2 c h{tn-tj)& h {f(U{)-f{u, h l{tj)))&S 
j=0 Jt i 

+ E f + A h^ / 2c h( t n-tj)& , h(f(u h ,l(tj)) - f{u h ,l{s)))ds 
j=0 Jt i 

+ E [ J A h~^ /2 ( C h(tn-tj)-C ll {t n -s))& , h f(u hi i(s))ds 
j=0 Jt J 

=: Ji + 72 + T 3 +/ 4 . 


Using the temporal regularity of the cosine operator, see (12. 1 II) . equation (12.131) . and the 
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assumptions (12 .8b (recall that g = 1 here), one gets 

E[||/i II 2 ] = I r \\A['- p \c h (t n - tj ) - C h (t n - s)) A^- 1)/2 ^||^o d, 
j=o J ‘j 2 

< c't [ ,J+1 ||A< - tj) c h (t n s)) AV-W&tA-W-W 

xA (P-i)/ 2 Q i/Y mds 
<C^min(2(^-l),l) for ^ e [1,2]. 

Next, using the convergence results from Theorem B.ll and the Lipschitz assumption on /, we 
observe that 


(E[||7 2 || 2 ]) 1/2 <'£ 1 f tj+i (E[\\A^~P )/ 2 C h (t n — tj)0 > h (f(U{) — f(ui h i(tj)))\\ 2 ]) 1 / 2 ds 
j= 0 J ‘j 

<c/£ 1 (E[||f//- M , il (t 7 )|| 2 ]) I/2 <Cfc min ^’ 1 ) for J8G[1,2]. 
j =o 


Similarly, using the assumptions on / given in (12.8b . and the regularity property of the finite 
element solution stated in Proposition ^. 31 one gets 

(E[||/ 3 || 2 ]) 1/2 <I / 0+1 -o)^(/(«/m(o))-/'(^, 1 (^)))ll 2 ]) 1 / 2 d.v 

j=0 J, J 

< ck min(P,l) for p g [1,2]. 

For the last term, J 4 , we obtain the estimate for /3 € [1,2] as follows 

(E[||7 4 || 2 ]) 1/2 < I I'*' (E[||A t P \ Ch (tn-tj)-C h (t n -s)) 

j=o J ‘J 


XA? 1)/2 ^ a / ! /(M/ i ,i(s))|| 2 ]) 1/2 dj- 




where we have used equation (12. 1 lb . the equivalence between the norms (12.14b . the assump¬ 
tions on the nonlinearity /, and the fact that the finite element solution u/ L 1 is bounded in the 

norm IMIm-i- 

Collecting all the above estimates and observing that 2(j3 — 1) < /3 for p £ [1,2], we 
finally get 


?|| 2 - ||«/ 1 ,2(f n )|| 2 ]| < for p e [1,2]. 

The second term in (14.5b can be estimated in a similar way as above. We have 

E[||Af t /ni 2 -iiAf M ,, 1 (f„)ii 2 ] <(E[||cr+ M , il (f n )ii 2 ^]) i / 2 (E[|| t /r- M/lil (f n )ii 2 2 ^]) 1 / 2 

< C(E[||Af-0 )/2 (f/f - M 1 ,*(r„))|| 2 ]) 1/2 

< cjt min(2(P-l),l) for p g [1,2]. 
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When estimating E[||a[~ — wi./ t (f H ))|| 2 ] we getthe same terms as J\ through /4 above, 

except that cosine is replaced by sine everywhere. Hence the same estimate holds. For the 
third and final term in ( 14.51 ). using the mean value theorem we get first 

M\V (t/f) - V K,i (*«)) IIlR®)] < CE[\\V(U?) - v (m /,.1 (f„))|| l 2 (@)] 

^ C\\V — Uf 1 \{t n )') Hl 2 (q,// 0 )• 

Recalling that f{u) = —V'(u), using Holder’s inequality, using the fact the numerical solu¬ 
tions are bounded in the mean-square sense, and the error bounds stated in Theorem l3. II we 
next estimate the following expression 

E[||V(t/f) — V(M/, i l(f„))|| il (g>)] < C\\v'— u h,i( t n))\\L 2 (n,H 0 ) 

< C (E[||t/f - u HA (t n ) \\l 2m i) 1/2 < cr(W). 

Putting all these estimates together we obtain equation ( 14.4b and the theorem is proven. □ 

5. Numerical experiments. This section illustrates numerically the main results of the 
paper. We first present the time integrators we will consider, then test their mean-square 
orders of convergence on various problems and finally illustrate their behaviours with respect 
to the trace formula from the previous section. 

5.1. Setting. The solution of our stochastic wave equation (11.11 ) will now be numerically 
approximated using the method of lines, i. e., with a linear finite element method in space and 
then with various time integrators (see below). Further, we will consider two kinds of noise: a 
space-time white noise with covariance operator Q = I and a correlated one with Q — A ' for 
some ,v > 0. We refer for example to 0 for a discussion on the approximation of the noise. 

We shall compare the stochastic trigonometric method (13.1b with the following classical 
numerical schemes for stochastic differential equations. When applied to the wave equation 
in the form ( 12 . 1 b . these numerical integrators are: 

1. The forward Euler-Maruyama scheme, see for example M or m, 

X n+l = X" + kAX” + kF(X n ) + G(X n )AW". 

2. The semi-implicit Euler-Maruyama scheme, see for example CD or m, 

X n +\ = x n + kAx n +1 + + G(X n )AW n . 

3. The backward Euler-Maruyama scheme, see for example Ifl4l or |fl 8 l , 

X n +\ = x n + kAx n +1 + kF (jn+ 1) + 

4. The semi-implicit Crank-Nicolson-Maruyama scheme, ini or |[ 20 l , 

X n +1 =x n + ^A(x n+l +X n ) + kF(X n ) + G(X n )AW n . 

Note that the backward Euler-Maruyama scheme, the semi-implicit Euler-Maruyama scheme, 
and the semi-implicit Crank-Nicolson-Maruyama scheme are implicit numerical integrators. 

All the numerical experiments were performed in Matlab using specially designed soft¬ 
ware and the random numbers were generated with the command randn. 
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FIGURE 5.1. The Anderson model: Spatial rates of convergence for the covariance operators Q = A 9 with 
5 = 0,1/2,1/3,1/4. The dotted lines are reference lines of slopes 1/3,2/3,10/18,1/2. M s = 2500 samples. 


5.2. Multiplicative noise. Let us first consider the one-dimensional hyperbolic Ander¬ 
son model ITIUOH: 

du(x,t) — Uxx(x,t)dt = u(x,t)dW(x,t) for (x,t) G (0,1) x (0,1), 
w(0,f) = n(l,f) = 0, t G (0,1), 

u(x, 0) = sin(27Dr), u(x, 0) = sin(37Tx), x G (0,1). 

This stochastic partial differential equation with multiplicative noise is now discretised in 
space by a linear finite element method with mesh size h. This leads to a system of stiff 
stochastic differential equations. The latter problem is then discretised in time by various 
integrators with time step k. 

Fiuu re [5~Tl i 1 lustrates the results on the spatial discretisation of the finite element method 
as stated in Theorem l3.ll The spatial mean-square errors at time 7) n( j = 1, 


Y^EfUu/Xx^end) -M(x,7;nd)|| 2 ] , 

are displayed for various values of the parameter h = 2~ l , t = 2,... ,9. The covariance op¬ 
erator is chosen as Q = A~ s for ,v = 0,1 /2,1 /3,1 /4. In the present situation, f(u) = 0 and 
g(u) = u satisfy the assumptions (12.81 ) with ft < s + J- This can be seen using the computa¬ 
tions done in Subsection 4.1 from MT3ll (with p = 2s and a = ^f-)- A clear dependence of 
the spatial convergence rates with respect to the covariance operator can be observed in this 
figure, in agreement with Theorem l3.ll Here, we simulate the exact solution u(x,t) with the 
numerical one using the stochastic trigonometric method (STM) (13.1b with a small time step 
£ ex act = 2~ 9 (in order to neglect the error from the discretisation in time) and /z eX act = 2~ 9 
for the mesh of the FEM. The expected values are approximated by computing averages over 
M s = 2500 samples. We computed the estimate for the largest standard errors, of all schemes, 
to be 0.0026. This shows that the error due to a Monte-Carlo approximation is negligible. 

We are now interested in the time discretisation of the above stochastic partial differ¬ 
ential equation with space-time white noise (Q = I and thus p <' 1 /2). We compute the 
temporal errors at time Tend = 0.5. In Figure [531 one can observe the rates of mean-square 
convergence of various time integrators. The expected rate of convergence ff(k l/2 ) of the 
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FIGURE 5.2. The Anderson model (space-time white noise): Temporal rates of convergence of the stochastic 
trigonometric method (STM), the semi-implicit Euler-Maruyama scheme (SEM) and the Crank-Nicolson-Maruyama 
scheme (CNM). The reference lines have slopes 1/4,1/3 and 1/2 . M s = 2500 samples. 


stochastic trigonometric method as stated in Theorem B.ll can be confirmed. Again, the exact 
solution is approximated by the stochastic trigonometric method with a very small time step 
^exact = 2“ u and uses /z exac t = 2“ 9 for the spatial discretisation. M s = 2500 samples are used 
for the approximation of the expected values. We computed the estimate for the largest stan¬ 
dard errors, of all schemes, to be 0.01. The numerical results for the forward and backward 
Euler-Maruyama schemes are not displayed since these numerical schemes would have to use 
very small time steps for such an /i eX act (see also Subsection 15.51 below). 

5.3. Semi-linear problem with additive space-time white noise. We next consider the 
sine-Gordon equation driven by additive space-time white noise (Q = I and thus P < 1/2) 

d«(x,f) — Uxx{x,t)dt = — sin(n(x,f))df + dW(x,f), (x,f) £ (0,1) x (0,0.5), 

«(0,f) = u(l,t) = 0, t £ (0,0.5), 

«(x, 0 ) = 0 , u(x, 0 ) = lji 3 j(x), x £ (0,1), 

where l/(x) denotes the indicator function for the interval I. 

Figure l573l disnlavs the rates of mean-square convergence at 7' enc j = 0.5 of various time in¬ 
tegrators. The expected temporal rate of convergence &(k 1 / 2 ) of the stochastic trigonometric 
method as stated in Theorem B.ll can be confirmed. Again, the exact solution is approxi¬ 
mated by the stochastic trigonometric method with a very small step size k exa ct = 2“ 11 and 
uses /i eX act = 2 9 for the spatial discretisation. M s = 2500 samples are used for the approx¬ 
imation of the expected values. We computed the estimate for the largest standard errors 
for all schemes to be 0.0027, showing that the error due to a Monte-Carlo approximation is 
negligible. 

5.4. Semi-linear equation with multiplicative noise. In this subsection, we consider 
the sine-Gordon equation driven by a multiplicative space-time white noise (Q = I and thus 

P < 1/2) 

du(x,t) — Uxx{x,t)dt = — sin(w(x,f))d? +u(x,t)dW(x,t), (x,t) £ (0,1) x (0,0.5), 

u(0,t) = u(l,t) = 0, f£ (0,0.5), 

u(x,0) = sin(27Tx), ii(x,0) = sin(37Tx), x £ (0,1). 
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FIGURE 5.3. The sine-Gordon equation (space-time white noise): Temporal rates of convergence of the 
stochastic trigonometric method (STM), the semi-implicit Euler-Maruyama scheme (SEM) and the Crank-Nicolson- 
Maruyama scheme (CNM). The dotted lines have slopes 1/4,1/3 and 1/2. M s = 2500 samples. 



FIGURE 5.4. The sine-Gordon equation with multiplicative space-time white noise: Temporal rates of con¬ 
vergence of the stochastic trigonometric method (STM), the semi-implicit Euler-Maruyama scheme (SEM) and the 
Crank-Nicolson-Maruyama scheme (CNM). The dotted lines have slopes 1/4,1/3 and 1/2. Af s = 2500 samples. 


Figure 15.41 displays the rates of mean-square convergence of various time integrators 
when applied to this semi-linear problem with multiplicative noise. The expected temporal 
rate of convergence of the stochastic trigonometric method as stated in Theorem l3.ll 

can be confirmed. One also observes a slower convergence rate for the other integrators. As 
before, a reference solution is computed by the stochastic trigonometric method with a very 
small step size /c eX act = 2“ u and uses /i eX act = 2~ 9 for the spatial discretisation. M s = 2500 
samples are used for the approximation of the expected values. We computed the estimate 
for the largest standard errors, of all schemes, to be 0.006. This shows that the error due to a 
Monte-Carlo approximation is negligible. 

5.5. Trace formula. We will now illustrate the trace formula from Section|4] To do this, 
we again consider the above sine-Gordon equation with additive noise and solve this problem 
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with a linear finite element method in space and in time we use the stochastic trigonometric 
method ( 13.11 ) with f{u) = — sin(«), g(u) = 1. Figure [53] (top) displays the expected value of 
the Hamiltonian along the numerical solutions of the above stochastic sine-Gordon equation 
where the covariance operator is given by Q = A 2 . In the present situation, the Lipschitz 
function /(«) = — sin(«) and the function g(u) = 1 satisfy the assumptions (12.81) with = 2 . 
This is seen using the fact that the eigenvalues of the Laplace operator with Dirichlet bound¬ 
ary condition satisfy A j ~ j 2 and the eigenvectors are given by {v / 2 sin(y' 7 tx)}j. The meshes 
are h = 0.1 and k = 0.01, the time interval is [0,5], and M s = 2500 samples are used for the 
approximation of the expected values. For this experiment, the largest standard errors for 
all the numerical schemes (except for the Euler-Maruyama scheme) is of the size of 0.002 
confirming that the Monte-Carlo errors are negligible. In this figure, one can observe the un¬ 
satisfactory behaviour of classical Euler-Maruyama-type methods. This is not a big surprise, 
since, already for stochastic ordinary differential equations, the growth rate of the expected 
energy along solutions given by these numerical solutions is incorrect H2lEl- The Crank- 
Nicolson-Maruyama scheme however seems to reproduce very well the linear drift in the ex¬ 
pected value of the Hamiltonian. Let us see what happens when one uses bigger time step and 
longer time interval. Figure [53] (bottom) displays the expected energies on the longer time 
interval [0,250] for the Crank-Nicolson-Maruyama and the stochastic trigonometric methods 
with a larger time step k = 0.1. The other parameters are the same as in the above numerical 
experiment, in particular, the Monte-Carlo error for the stochastic trigonometric method is 
negligible (error of size 0.01). On this long-time interval, excellent behaviour of the stochas¬ 
tic trigonometric method ( 13.11 ) is still observed although this does not follow from the result 
presented in Theorem l4.5l 

6 . Appendix. In order to improve the readability of the paper, we give some details for 
the proofs of the results given in Section [2] The proofs of ( 12.91 ) and (12.10b can be found in 
Corollary 4.2 in fl 6 l . They are obtained by interpolation between the results for the endpoints 
of the parameter values. These in turn are well-known estimates for the finite element approx¬ 
imation of the homogeneous wave equation: «/, + A/,m/, = 0, t > 0; m/,(0) = &i,uo, «/,(0) = 
&hV o- 

For example, for 7 = 1 we have, by a standard stability estimate, 

\\m)Xo\\ < C (II^Koll + ll^voll-t.fc + IM + 11 VO 11 — 1 ) 

<C(||«o||i + ||vo||-i)<C/,°|||Xo||| 1 , 

since /A), is not bounded with respect to the //°-norm. For y = 3, we have 

II W*o|| < C(1 +t)h 2 (IMb + llvollz) < C( 1 +f)fi 2 |po||| 3 , 

cf. the estimation of f), in the proof of Corollary 4.2 in m. Interpolation between these two 
cases completes the proof of the first bound in (12.9b . Note that the required initial regularity 
is one order higher than the order of convergence. This is typical of the finite element method 
for the wave equation. Another choice of projector, m/,(0) = S?hUa, would give a slightly 
better result for low initial regularity here, but a worse result for Sf/,. 
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